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A new kind of dynamo utilizing flowing laboratory plas- 
mas has been identified. Conversion of plasma kinetic energy 
to magnetic energy is verified numerically by kinematic dy- 
namo simulations for magnetic Reynolds numbers above 2f0. 
As opposed to intrinsically-turbulent liquid-sodium dynamos, 
The proposed plasma dynamos correspond to laminar flow 
topology. Modest plasma parameters, f-20 eV temperatures, 
10 19 -10 20 m -3 densities in 0.3-1.0 m scale-lengths driven by 
velocities on the order of the Alfven Critical Ionization Veloc- 
ity (CIV), self-consistently satisfy the conditions needed for 
the magnetic field amplication. Growth rates for the plasma 
dynamos are obtained numerically with different geometry 
and magnetic Reynolds numbers. Magnetic-field-free coaxial 
plasma guns can be used to sustain the plasma flow and the 
dynamo. 

PACS numbers: 52.30.-q; 47.65. +a; 52.65. Kj; 52.72,+v 



Dynamo action is believed to be the fundamental 
mechanism that creates the magnetic field commonly ob- 
served from planetary to galactic scales ||. Conversion 
of kinetic (flow) energy into magnetic energy is the best 
known explanation for observed magnetic fields in the 
universe that are persistent against resistive diffusion. A 
direct experimental proof of the conversion of plasma flow 
energy to magnetic energy has not been accomplished to 
date. Another motivation of this letter is to show dy- 
namo experiments can be performed in plasmas, as an 
alternative to liquid sodium or other conducting liquids 
0. The plasma dynamos discussed here are funadmen- 
tally different from what has generally been known for 
spheromaks |||| and reversed field pinches [0J|], when 
the term 'dynamo' has traditionally been invoked to ex- 
plain conversion of one type of magnetic flux into another 
(toroidal flux into poloidal flux in the spheromak case). 
Because in these previous cases, the total magnetic field 
energy does not increase. In the present case, the plasma 
flow energy is converted into magnetic energy, and mag- 
netic field is amplified accordingly. 

Existing and proposed laboratory dynamo experiments 
|^| |]T| have used liquid sodium as conducting medium. 
Liquid sodium has low resistivity and viscosity fl2]| for 
conduction of electricity and fluid flow, as well as a 
low melting temperature which eases the creation of the 
liquid state for the conductor. A constraint on liquid 
sodium dynamos is intrinsically turbulent flow due to the 
small ratio of viscosity to resistivity jl^] . The turbulent 
flow makes it difficult to make detailed comparisons be- 
tween experiments and theories. Also, resistivity of the 
liquid sodium is only variable within a factor of two p0[ , 



which implies the magnetic Reynolds number, Re m , may 
be varied only within a small range, primarily by varia- 
tion of the experimental dimensions. Even with the high 
conductivity of the liquid sodium, existing technologies 
can only move liquid sodium at velocities no more than 
20 m/sec jl^]. Therefore ~1 meter in size is necessary 
to have a dynamo excitation for a liquid sodium exper- 
iment, and Re m has been limited to values below a few 
hundred 

Dynamo action is described by the induction equation 
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where U is the fluid flow velocity, and rj is the resistivity. 
Without the flow U, the induction equation is a diffusion 
equation for magnetic field B with the characteristic dif- 
fusion time determined by the dimension of the field and 
resistivity. Only when U does not vanish may B grow 
by transfcring the flow energy. The magnetic Reynolds 
number Re m = fioLUo/rj measures the relative ampli- 
tude of the flow drive to diffusion in Eq. (|l|), where L, 
Uq and rj are characteristic values for the dimension, ve- 
locity field, and resistivity. Theories |l2|,[l3| predict that 
only when the magnetic Reynolds number Re m exceeds 
certain threshold values (from 20 to 200, depending on 
the geometry), is dynamo action possible. 

Re m greater than 200 appears easily achievable in flow- 
ing laboratory plasmas by this analysis. Thus an ex- 
periment is possible to create magnetic field energy that 
grows by transformation of the plasma flow energy alone. 
Such a new plasma dynamo can cover a wide range of 
magnetic Reynolds numbers both below and above the 
threshold by controlling the plasma flow velocity and the 
plasma resistivity. In the entire plasma system, one may 
encounter various length scales L, velocities U, as well 
as uncertain resistivity (the plasma resistivity is often 
"anomalous"). Average, "global" values characteristic of 
the entire flow appear reasonable to use for comparison 
to theoretical threshold values of Re m . In the plasma 
regime considered here where the flow is laminar (see 
below), the scale length L is close to the dimension of 
the boundary. The characteristic velocity Uq is repre- 
sentative of a Beltrami-fiow profile pi} , which is a flow 
satisfying V x U = AgU with constant parameter A#. 
Spitzer resistivity provides a realistic estimate since the 
magnetic field is weak during the initial growth of the 
magnetic field, and the long plasma lifetimes imply local 
thermal equilibrium. In addition, energy balance is dom- 
inated by atomic physics (radiation by the electrons and 
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boundary losses by the ions) and classical thermal con- 
duction ||15|j , and the plasma temperatures will be close 
to being in equilibrium and in the 1-20 eV range, even 
with substantial high-Z impurities. 

One key feature of these laboratory plasma dynamos 
is the high speed plasma flows (see Fig. 1) which will be 
on the order of the critical ionization velocity (CIV) p6[ . 
CIV, denoted as U c , is defined as U c = y/2Ei/Mi, where 
Ei is the ionization energy and Mi the ion mass. A veloc- 
ity of this magnitude is routinely achieved using plasma 
accelerators or thrusters |17j]. For hydrogen gas, U c is 51 
km/sec. For ease of estimation, the electron temperature 
determining the resistivity can be assumed to be in equi- 
librium with the ion temperature which is itself set to the 
"flow equivalent temperature" Tj = (MiU 2 )/2. With this 
simplifying assumption, Re m scales as Re m — 15T? [eV] 
L [m]. For L about 0.5 m and Ti ranging from 1 to 
20 eV, the Re m varies from 7.5 to 3000. The threshold 
Re m — 200 corresponds to T,; = 5.2 eV, or a flow velocity 
of 31 km/sec for hydrogen, and hence in the expected 
range of CIV velocities. Because of the high plasma flow 
velocity, Re rn in a plasma dynamo can be well above 
the threshold magnetic Reynolds number for magnetic 
field to grow, even considering that the resistivity of the 
plasma is much greater than that of liquid sodium. In 
Fig. 1, these plasma dynamos are compared with exist- 
ing and proposed sodium dynamos, which are the only 
known laboratory dynamos so far. 

Another unique feature of the new plasma dynamo is 
laminar flow. Flow topology is determined by the ki- 
netic Reynolds number, Re, which characterizes the rel- 
ative amplitude of the fluid convection to viscous dis- 
sipation in the Navier-Stokes equation. The critical 
Reynolds number for onset of turbulence Re is more 
than a few thousand, and fluid motion is laminar for 
Re below this threshold Jig] , Plasmas discussed here 
can be approximated as fluids because the mean free 
path for particle motion is much less than the system 
dimension. The kinetic Reynolds number Re then scales 
with the device dimension and plasma ion temperature 
as Re = 4.0 x 10~ 16 L [m] n e [m- 3 }/T 2 [eV]. An ex- 
perimental dimension L of 0.5 m is adequate to meet 
the conditions of a plasma dynamo, with density 10 19 to 
10 20 m~ 3 and temperature 1-20 eV. This corresponds to 
the kinetic Reynolds number Re varying from 5 to 20000, 
which imply that plasma flow will be mostly laminar. 

Plasma flow can be created using coaxial plasma guns, 
or Marshall guns |p|Jl9|. For the gun plasma momen- 
tum to be effectively transferred to the bulk plasma in a 
chamber, the mean free-path, A m / P , of the gun plasma 
ions should be less than the size of the dynamo cham- 
ber, that is, X m fp < L is required. This is equivalent to 
Re > 2.5, and hence the kinetic Reynolds number can- 
not be too small. The magnetic Prandtl number (Pr m ) 
is defined as the ratio Pr m = Re m /Re and scales as 
Tfjn e with values in the range of 3.7xl0~ 4 to 600 for 
the plasma parameters above. In comparison, the Pr m 



for liquid sodium experiments is ~ 10 -5 . The plasma 
Prandtl number can be very similar to that of the inte- 
rior of the sun or the galactic plasma dynamos, albeit at 
very much lower Reynolds numbers. 

The gun-plasma momentum transfer is also affected 
by ion- neutral collisions, which include ion-neutral mo- 
mentum transfer and ion-neutral charge exchange. The 
ion-neutral momentum transfer cross section is 100 times 
less than that of ion- ion collision |2(J. Therefore, as long 
as the ionization fraction is greater than 1%, ion-neutral 
momentum transfer can be neglected. The charge ex- 
change cross section is about half of that of ion-ion mo- 
mentum transfer. Therefore ion-ion momentum transfer 
is the dominant process for gun-plasma momentum trans- 
fer. Generation of toroidally rotating plasmas is most ef- 
ficient when the ionization fraction is greater than 30%. 

A possible experimental configurationtion is shown in 
Fig. 2, when one 'main gun' (Labeled as MG) is used to- 
gether with three 'toroidal guns' (TG, not shown) to pro- 
duce toroidal and poloidal plasma flows. Kinetic energy 
density for a 13 eV and 10 19 m~ 3 plasma is equivalent to 
a magnetic field strength of 53 gauss at the same energy 
density. Theoretical studies indicate that the fraction 
of the plasma flow kinetic energy that can be converted 
into magnetic field depends on the ratio of the magnetic 
Reynolds number to the kinetic Reynolds number [El]. 
As much as 20% of the flow energy was seen to con- 
vert into magnetic field numerically. If this is the case, 
magnetic field of up to 24 gauss may be generated when 
the back-reaction of the magnetic field on plasma flow is 
taken into account. 

Numerical results confirm the existence of laminar 
plasma dynamos. The plasma flow is approximated by 
a Beltrami flow, V x U = AgU, which corresponds to 
a state of maximum kinetic helicity [J dV(XJ ■ V x U)[ 
for a given total kinetic energy and a given boundary 
condition. We can use the Beltrami flow approximation 
because the flow is laminar here, large coherent flow pat- 
tern can be established. The Beltrami flow is also an 
eigenmode for the cylindrical boundary condition con- 
sidered. Theoretical results indicate that large kinetic 
helicity content can lead to large dynamo growth rate. 
The solution of axisymmetric Beltrami flow inside the 
cylinder < r < R and — L/2 < z < L/2 may be written 
in terms of a poloidal flux function \I/(r, z) 

(V \ 7TZ 
iii^Jcos— , (2) 

„._!£, (3) 

r oz r or r 

where J\ (x) is a first order Bessel function in its standard 
notation, and jn is the first root of J\{x) — 0, \ 2 B = 
j\\IR 2 + tt 2 /L 2 . The velocity components U r , U z , and 
U,), have been normalized to a characteristic velocity Uq. 

The kinematic dynamo problem for Beltrami flows sat- 
isfying Eqs. (||) and (^) is implemented numerically as 
follows. The flow field is taken as given and the back 
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reaction of the growing magnetic field on the flow is ne- 
glected. This approximation is justified for the initial 
stage of the exponential dynamo growth, when the mag- 
netic field is weak and does not influence the flow. In- 
stead of solving the induction equation Eq. ([!]) for B 
directly, potentials A and ip are introduced. Using the 
gauge condition ip — U • A+ (r)/ fj,o)V ■ A = |22|], one can 
derive the following equation of evolution for the vector 
potential A, 

dA l = A dU k lTk dA> | r, 8 2 A* 
dt dx % dx k no dx k dxk ' 

where the resistivity rj is assumed to be constant through- 
out the cylinder and the coordinate notations refer to a 
Cartesian coordinate system x l . A 3D kinematic dynamo 
code is used to solve Eq. (Q). Then, the magnetic field 
can be obtained at any time by taking the curl of A. The 
code is written for cylindrical coordinates, and it uses an 
explicit scheme with central spatial differencing in the 
advection term and standard nine points stencil for the 
diffusion term. Since the conductivity of metallic walls 
is much higher than the conductivity of the plasma, a 
perfectly conducting boundary is a good approximation. 
All boundaries of the cylinder (r = R, z = —L/2, and 
z = L/2) are assumed to be perfect conductors. Then, 
the boundary conditions for A at perfectly conducting 
boundaries compatible with the gauge can be chosen as 
follows: both the components of A parallel to the bound- 
ary and the divergence of A at the boundaries are zero. 
This gives three boundary conditions for three compo- 
nents of vector potential. Eq. (||) has a unique solution. 
A detailed description of the code, of the gauge choices, 
and of the implementation of the boundary conditions 
can be found in Ref. [p3| . 

In the case of axisymmctric flow, the nonaxisymmet- 
ric modes of the magnetic field, which are proportional 
to exp(m<^ — iu)t) with different azimuthal wavenumbers 
n (also known as toroidal mode numbers), are decou- 
pled from each other and are eigenmodes of Eq. (|l|) or 
Eq. ([|). The 3D kinematic code picks up the fastest 
growing mode of the dynamo, which turns out to be the 
n = 1 mode. Since the structure of each eignemode is 
two dimensional in (r, z) coordinates, a 2D code involv- 
ing only the azimuthal component of the magnetic field 
saves much time for simulations with varying boundary 
parameters. Such a 2D code for the vector potential evo- 
lution has also been written and used to calculate the 
dependencies of the growth rates of the n = 1 mode 
on the boundary dimensions. The growth rates and the 
structure of the n = 1 modes obtained using the 3D 
dode agree remarkably well with that using the 2D code. 
The modes have an oscillatory nature, i.e. the oscilla- 
tion frequency w — ui' + ij(U max /L) has real and imag- 
inary parts, where both the rotation frequency to' and 
the dimensionless growth rate 7 are real. U max /L is 
about 100 kHz. U max is the maximum absolute veloc- 
ity inside the cylinder. Exponentially growing (or decay- 
ing) magnetic fields also rotate with certain frequency, 



which comparable to the frequency of the fluid rotation, 
J ~ 100 kHz. 

Dependence of the growth rate 7 on the magnetic 
Reynolds number Re m and on the aspect ratio R/L of 
the cylinder is explored. Fig. 3a shows the dependence 
of 7 on Re m for a fixed aspect ratio R/L = 1. Nega- 
tive values of 7 correspond to decaying magnetic field. 
For Re m smaller than the threshold value of 210, the dy- 
namo will not grow. With increasing Re m , the growth 
rate first increases up to the maximum value ~ 0.022 (a 
growth time of 0.45 msec) at Re m = 500, and it subse- 
quently decreases to smaller values but remains positive. 
Such behavior is typical for slow dynamos (for example, 
see Ref. ||). Thus, these laminar plasma dynamos are 
slow dynamos, for which 7 asymptotically goes to zero 
for Re m — > 00. The laminar and regular (as opposed 
to chaotic) fluid motion stretches the magnetic field lin- 
early in time. Exponential growth is only possible due 
to small diffusivity of the magnetic field, which recreates 
the component of the magnetic field perpendicular to the 
flow. When the diffusion of the magnetic field becomes 
very small (R m — > 00), the growth rate 7 also decreases 
to 0. Fig. 3b shows the dependence of 7 on the aspect 
ratio R/L for a fixed Re m = 500. Dynamo does not ex- 
ist for very long or very short cylinders. It is necessary 
to have comparable R and L for efficient excitation of 
the dynamo. In particular, at an experimentally feasible 
value of Re m = 500, the maximum growth rate 7 = 0.055 
is achieved at R/L = 0.6. These results suggest that a 
somewhat elongated cylindrical vessel will be the best 
for the excitation of dynamo. Pulsed plasma flow of sev- 
eral msec long should be sufficient to excite the laminar 
dynamos. 

In conclusion, a theoretical study has shown dynamo 
can be excited to convert plasma kinetic energy into mag- 
netic field energy in a laboratory environment. These 
plasma dynamos are laminar because of the low ki- 
netic Reynolds number. Modest plasma parameters self- 
consistently satisfy the conditions needed for the dy- 
namo. Numerical calculation yields a threshold magnetic 
Reynolds number of 210 for exponential growth of the 
laminar plasma dynamos in a cylindrical boundary with 
Beltrami flows. These results indicate that a new type of 
laboratory dynamo experiments is possible. 
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FIG. 1. Flow velocity vs experimental size for various 
sodium dynamos and for laminar plasma dynamos (Shaded 
region). Specific region of the laminar plasma dynamos may 
vary with gas species. 
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FIG. 2. A possible laminar plasma dynamo configuration 
using coaxial plasma guns. 
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FIG. 3. The dimensionless growth rate of the toroidal 
mode n = 1 of a laminar plasma dynamo (7) as a func- 
tion of the magnetic Reynolds number, Re m , for aspect ratio 
R/L = 1. 
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Dependence of the dimensionless growth rate of 
the toroidal mode n = 1 of a laminar plasma dynamo (7) on 
the aspect ratio, R/L, for fixed magnetic Reynolds number 
Re m = 500. 
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